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way to study the phase structure. We estimate the location of the critical point at finite chemical 
potential in the heavy quark region. 
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1. Introduction 

In our previous study on the QCD phase structure at zero chemical potential in the heavy quark 
region we have identified the phase structure through the shape of an effective potential Veff (P) 
defined by the probability distribution function of the plaquette P. Since the plaquette is equal to 
the gauge action divided by the gauge coupling parameter j8 , we can precisely evaluate Vgff in a 
wide range of P by combining data at different j8 's via the reweighting technique. 

For the same method to work at non-zero chemical potential and/or small quark masses, 
reweightings in the direction of the quark mass and the chemical potential are also mandatory. 
Therefore, it is useful to introduce new Veif which depends on variables other than the plaquette. A 
study in a light quark mass region by this method is reported in |^. 

In this report, we consider the phase structure at non-zero chemical potential in the heavy 
quark region. At the leading order of the hopping parameter expansion, the influences of dynamical 
quarks are given in terms of the Polyakov loop Q.. Using a new Veff which depends on both P and 
= ReQ., we study the QCD phase structure at zero and non-zero chemical potential. We adopt 
the cumulant expansion method ||3|, 0], in order to control the sign problem due to the fluctuation 
of the complex phase of the quark determinant at finite chemical potential. The method with new 
Veff works well to study the fate of the first order transitions. We evaluate the location of the critical 
point where the first order deconfining transition in the heavy quark limit turns into a crossover, at 
zero and non-zero values of the chemical potential. 

2. The method 



In this paper, we employ the standard plaquette gauge action and unimproved Wilson quark 
action, given by 



-6N,it,pP, 

f=l I n n 1^=1 
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where P=-f^ £ Re Tr Un,^iU„^^yUl^~^^Ul„ is the plaquette, jS = is the gauge cou 



n,/i<v 



pUng parameter, Kf {jXf) is the hopping parameter (chemical potential) for the /-th flavor, and 
^site = X is the lattice volume. In this study, we mainly consider the degenerate case Kf = K 
and Hf = n for f = I, - ■ ■ ,Ni. Note that M,,^ we study does not depend on j3. 
The probability distribution function for P and Q.r = R&O. is defined by 



wiP,aR;l5,K,^) = J m5{P[U]-P)5{D.K[U]-aR)[detM{K,^)f'e^^-''P'', (2.3) 
D.R = Ren, D. 
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We now define the effective potential as 
Applying the reweighting technique, we find 



w{P,^r;PAO) 



(2.4) 



(2.5) 



where Vo{P,D.r;P) = -lnw{P,D.R;p,0,0) = Veff (PjHr; j8,0,0) is the effective potential in the 
heavy quark limit. The second term in the r.h.s. is the reweighting factor, which is actually in- 
dependent of p as 



w{p,aR-p,K,^) l!^u5iP-p)5{CiR-aR) 



detM{K,ll) 
detM(0,0) 



w{P,Q.r;P, 0,0) 



jm5{p-p)d{aR-aR) 



detM{K,iJ.) 
detM(0,0) 



Since the j8 -dependence is inherent in Vq only, reweighting in j8 is simply given by 
V,s{P,aR;l5,K,^) = V,siP,QR-l5o, k,^) - 6N,it,il5 - po) P 



(2.6) 



(2.7) 



At the lowest order of the hopping parameter expansion, the ratio of the quark determinants in 
the r.h.s. of (12.61) is evaluated as 



detM(K-,jU) 
detM(0,0) 



: exp 



2SSN,ir^K^P + 3 X 2^'+^N^K^' jcosh (^^) + /sinh (^^) Hi} 



(2.8) 



where flj = ImQ. is the imaginary part of the Polyakov loop. The first term proportional to P can be 
absorbed into the gauge action by a shift j8 — j8* = j8 +48A^fK"^. In the evaluation of the effective 
potential at non-zero chemical potential, the most difficult part is to evaluate the expectation value 
of the complex phase at fixed values of P and Q.r, denoted by 



e'" ) with e = 3 X l^'+^NfN^Xai, X = K^' sinh [l^/T). 



(2.9) 



3. Results at 11 = 



Let us first discuss the case of = 0, in which the imaginary phase term in ( p^q ) is absent. 
The reweighting formulae for the derivatives of Veff simply become 



^{P,a^;l5,K) = ^(P,a^-p,)-6N,i,,{l5+4mK'-l5o), 



da 



R 



-{p,^r;k) 



do. 



(3.1) 
(3.2) 



R 



where the argument j8 in dVsff/dQ.R and dVo/dO-R is omitted in (3.2) since they are independent 



of j8 due to (p77|). Note that, besides overall constants [the last terms in (3.1) and (3.2)], the de- 
pendences of dVsn/dP and dVefi/dQ.R on P and D.r are independent of j8 and K [the first terms in 
(0) and (0)]. 
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In Figs. [T] and ^ we plot dVo/dP and dVo/dO.^ as functions of (f,nR). Data are taken from 
the pure gauge configurations generated on a 24^ x 4 lattice in [|l|]. To evaluate w{P, Hr; j8 , 0, 0), we 
replace the delta function by the Gaussian function: d{x) ^ exp [— (a;/A)^] where A = 0.0005 
for P and 0.01 for Hr. Using we combine data at five jS's in the range jS = 5.68-5.70 to 
improve the statistics in a wide range of P and Hr. We also average D. over the Z(3) rotation, 
and remove statistically poor data points with w < 10^^^^ x JwdPdQ.R. We then calculate the 
derivatives by the difference between x — e/2 and x + e/2, where £ = 0.00025 for P and 0.005 for 
Hr. These parameters are determined by consulting the statistical tolerance of final results. 

To study the phase structure, it is useful to follow the curves dV^s/dP = and dV^fi/dO.^ = 0. 



From ( |3.1| ) and (p^, these curves at different (j3,K") corresponds to different contour curves of 
dVa/dP and dVo/d^R, as plotted in Fig. |. 

When the curves dV^n/dP = and 5Veff/^^R = cross at only one point, we have just one 
minimum of Veff- In this case, the vacuum changes continuously by a small shift of j8 and K, 
indicating no first order transitions around this (j3,K"). On the other hand, when we have three 
intersection points (see the contour curves dVo/dP = —1000 and dVo/dQ.R = 50 as well as the 
thick stars in Fig. ^), we have two minima and one saddle point, implying the existence of the first 
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Figure 4: Critical point in the phase- Figure 5: (S^) and (0"^)^, at A = 2.0 x 10"'' around the 

quenched approximation for N{ = 2. critical point {P, Q.r) « (0.546, 0.06) for A^f = 2. 



order transition. In particular, from the merger of three intersection points, we can determine the 
critical point where the first order transition line terminates. 

From Fig. ^, we find that three intersection points appear due to the S-shape of the curve 
dVeff/dQ.^ = at small K. The S-shape becomes weaker with increasing K, and eventually the 
three intersection points merge to one intersection point at the critical point (e.g. the intersection 
of the contour curves dVo/dP = —4000 and dVo/dQ.^ = 150 shown by a thin star in Fig. ^ seems 
to be beyond the critical point). By consulting the contour curves, we find that dVo/dP —3150 
and dVo/dD.R 120 (illustrated by bold curves in Fig. |3|) seem to correspond to the critical point. 
This gives us a preliminary estimation K^p 0.0690(7) and j8*p « 5.6805(2), corresponding to 
j3cp ~ 5.6783(2) for the case of N{ = 2, where the errors are estimated from the error of the contour 
curves due to the statistical errors of potential derivatives. 

Our previous study using the effective potential for P gives K^p = 0.0658(3)(lj[ ) and j8*p = 
5.6836(1) (5) (j3cp = 5.68 19(1) (5)) for A/^f = 2, where the second errors are systematic ones due to 
the method-dependence of the results Our new estimation deviates slightly from the previous 
one. We are currently testing a refinement of the method to extract smoother contour curves. In 
the remaining part of this report, however, we instead proceed to the case of / 0, by taking an 
advantage of the relatively simple reweighting formulae for jj.. 



4. Results at /i 7^ 



At jti / 0, we need to evaluate the complex phase factor /pQ^ defined by ([2.9|). When 6 



fluctuates a lot at large jU., it becomes difficult to estimate )pQ^ reliablly (the sign problem). 

4.1 The case of phase-quenched finite density QCD 

We first consider the case of phase-quenched finite density QCD, in which the complex phase 
term is removed in the quark determinant ratio ( |2.8[ ). In two-flavor QCD, this corresponds to the 
case of the isospin chemical potential, /i„ = —jx^ = After absorbing the first term in the r.h.s. of 



(2.8) into the gauge action by j8 — )• j3*, we find that the effect of fX is just to modify the hopping 



parameter as k: — )• k cosh^^^' (n/T) in the theory at /i = 0. Therefore, to the lowest order of the 
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hopping parameter expansion, the critical point at finite jj. is given by 

K^p(At) = Kip(0)/cosh'/^' (At/r), (4.1) 

where JQ;p(0) is the critical point at vanishing chemical potential. See Fig. Note that, with 
increasing /i, the critical point approaches towards K = where the hopping parameter expansion 
is reliable. 

4.2 Cumulant expansion 

We now turn on the complex phase. In order to incorporate the effects of the complex phase, 
we adopt the cumulant expansion method M, H] : 



= exp 



(4.2) 



where d is given by and (0)^. = (0), (e^)^ = (^2^ _ {ef, (03^^ = ^^3^ _ 3 ^^2^ ^ 
2 (0)^, etc., with the averages taken at fixed P and IIr. Since the symmetry of QCD under 6^—6 



implies {d^"^^)c = 0, the r.h.s. of (4^) is manifestly real and positive. This means that, if the 
cumulant expansion converges, the sign problem is solved. 

The distribution of the complex phase in the quark determinant has been shown to be well 
approximated by a Gaussian distribution up to moderate values of jj. when the complex phase 6 
is appropriately defined ^, This means that the leading order term {d^)^ dominates in the 
cumulant expansion, and hence the expansion converges. In the heavy quark region, we study 
(02«^>^ = (3 X 2^'+^N{N^Xf" {^i")^. with A = k:^' sinh(At/r). Our results for (e^)^ and the next- 
leading order (0^)^ are plotted in Fig. ^ for A = 2 x 10^^. As will be discussed in the next 
subsection, this A corresponds to the critical point in the large /i limit. Note that the scale for 
{0^)^. is much magnified in this figure. We find (6^)^ ^ (^'*)c' accordance with the Gaussian 



dominance. In the foUowings, we take the leading order approximation of ( [4.2| ) assuming the 
Gaussian dominance. 

4.3 Critical curve at / 

At the leading order of the cumulant expansion, derivatives of Vgff read 

^^^^ f.N (R^Am.^ R) , (3x2^'+^M-A^3^)^^(^?). 

^ = ^ - 6A^site [P +4m{K - l3o) + 2 ^^-^^ 
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When the last term in ( |4.4D with the factor d {^j)^ /^^r modifies the S-shaped contours of dVo / 
shown in Fig. ^ the critical point shifts from that of the phase-quenched case given by (11). Our 
results ofd{D.j)JdP and d{Q.j)JdaK are shown in Fig. n. 

We note that d(Q^^ ^JdP is quite flat and small around the critical point, and thus just causes a 
shift of jScp according to ( [4.3[ ). We also find that the effect of d {^)^ / dQ.R is quite small around the 
critical point: When we disregard the last term in ( |4.4| ), the critical point locates at k:'^' cosh(/x/r) = 
jQ;p(0)^' ^2x 10"^ as discussed in Sect, ^j], i.e. X^2x 10"^ x tanh(/x/r) along the critical curve. 
Using these values and consulting the right panel of Fig. ^, we find that the contribution of the last 
term in ( p^ ) is at most about 3% of the second term around the critical point even in the large /i 
Umit where tanh(/x/r) 1. Therefore, KcpilJ.) is indistinguishable from K;fp(jU) shown in Fig. 0, 
within the current statistical errors. 



5. Conclusion 

We have studied the phase structure of QCD at non-zero chemical potential /i in the heavy 
quark region, using the effective potential defined by the probability distribution function of the 
plaquette P and the real part of the Polyakov loop Q.^. The reweighting technique enables us 
to obtain the effective potential at 7^ in a wide range of P and Hr. Adopting the cumulant 
expansion method to calculate the effects of the complex phase, we have shown that the derivatives 
of the effective potential provide us with an intuitive and powerful way to investigate the fate of 
first order phase transitions. 

We find that the critical point where the first order deconfining transition in the heavy quark 
limit terminates locates quite close to that in the phase-quenched case up to large values of jj.. 
The smallness of the effects of the complex phase around the critical point is due to the fact that 
Kcp becomes rapidly small as jj. is increased. Therefore, a careful examination of the effects of 
the complex phase is required off the critical region and at lighter quark masses. An attempt to 
study finite density QCD at light quark masses by combining phase-quenched simulations and the 
reweighting technique is reported in [Q]. 
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